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Abstract  Burlington  VT  Background 


This  paper  discusses  the  aerodynamic 
predictions,  tests,  and  analyses  of  a  slender  fin 
stabilized  missile  configuration  at  Mach  numbers 
from  1.0  to  5.0.  Prediction  techniques  consisted 
of  both  empirical  and  analytical  methods, 
including  state-of-the-art  computational  fluid 
dynamic  codes.  Free  flight  tests  were  also 
conducted  on  subscale  models  to  obtain  an 
aerodynamic  baseline  to  which  the  various 
predictions  could  be  compared.  This  paper 
summarizes  these  results,  attempts  to  identify  the 
physics  governing  the  flight  trajectories,  and 
proposes  new  ideas  for  further  research  in  this 
technology  area. 

Nomenclature 

p 

A  =  reference  area,  d  / 4 

Clp  =  roll  damping  derivati”e 

Cma  *  Pitchi-ng  moment  derivative 

Cmq  =  pitch  damping  derivative 

Cu  =  normal  force  derivative 

CXo  =  zero  yaw  axial  force  coefficient 

d  =  projectile  body  diameter  (reference 

length) 

1  =  projectile  length 

M  =  Mach  number 

p  =  3pin  rate 

Rj^  =  Reynolds  number  based  on  length 
v  =  velocity 


Fin  stabilized  missile  configurations  are 
known  to  exhibit  decreasing  static  stability  as  a 
function  of  increasing  Mach  number.  Figure  1 
shows  a  typical  static  stability  plot  where  the 
missile's  restoring  moment  continuously  decreases, 
eventually  leading  to  static  instability  (l.e.  the 
pitching  moment  slope  becomes  positive).  Because 
of  the  obvious  advantages  associated  with 
increasing  the  3peed  of  missile  configurations, 
the  Air  Force  Armament  Laboratory  (AFATL) 
initiated  an  in-house  program  to  investigate 
potential  methods  to  extend  future  flight  regimes. 
However,  before  this  could  be  accomplished  It  was 
believed  that  a  substantially  more  detailed 
understanding  of  the  flow  physics  creating  this 
situation  must  be  obtained.  Therefore,  the  first 
phase  of  this  program  was  to  predict, 
substantiate,  and  understand  the  static  stability 
characteristics  of  a  generic  missile 
configuration.  This  paper  summarizes  the  current 
research  which  attempts  to  fulfill  those  goals.  A 
corollary  to  this  research  is  to  begin  identifying 
key  missile  geometry  parameters  (l.e.  fin 
shape/size)  which  will  effectively  increase  the 
stability  envelope. 


X  ;  projectile  longitudinal  axis 
Xcg/1  =  center  of  gravity  location 
Xcp/1  —  center  of  pressure  location 
&  =  effective  angle  of  attack  squared 
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Fig.  1  Typical  Cma  versus  Mach  Number  plot. 

Test  Configurations 

In  determining  a  configuration  on  which  to 
conduct  this  research,  it  was  decided  to  select  a 
configuration  which  would  most  dramatically 
exhibit  the  static  stability  limitations 
previously  described.  Therefore,  a  generic,  high 
fineness  ratio,  tail-controlled  missile 
configuration  was  developed  as  the  baseline  model 
(Figure  2).  The  model  consisted  of  an  ogive  nose, 
cylindrical  body,  with  one  four-fin  group  located 
near  the  extreme  aft  end  of  the  projectile.  Two 
separate  fin  configurations  were  fabricated  for 


Test  Facilities 


the  free  flight  tests,  the  difference  being  that 
the  delta  fin  planform  area  of  one  configuration 
was  exactly  twice  the  area  of  the  other 
(maintaining  the  fin  aspect  ratio).  All  fins  were 
equally  spaced  at  90  degrees.  The  fin  thickness 
(made  of  .015  inch  shim  stock)  was  small  compared 
to  the  body  diameter.  The  fabricated  models  were 
approximately  5  percent  scale  compared  to  current 
air-to-air  missile  configurations  and  had  a 
fineness  ratio  (length  to  diameter)  of  29  which  is 
representative  of  today's  air-to-air  missiles.  In 
addition  to  the  two  fin  designs,  models  consisting 
of  two  separate  center  of  gravity  locations 
(approximately  93  and  52  percent  length  from  the 
nose)  were  tested  for  each  design.  All  models 
were  constructed  in-house  and  passed  a  quality 
control  inspection  prior  to  testing. 


-Ml.Mh-  1 


Fig.  2  Sketch  of  configurations. 

The  model/3abot  pickages  were  designed  to  be 
fired  from  a  smoothbore  30mm  barrel.  A  single 
stage  powder  gun  was  used  to  obtain  flight  data  up 
to  Mach  3.5  and  a  two-stage  light  gas  gun1  was 
utilized  to  achieve  the  flight  data  up  to  Mach 
5.0.  Figure  3  is  a  photograph  of  two  model/sabot 
packages,  showing  both  the  large  and  small  fin 
models.  The  sabot  serves  several  Important 
functions.  It  protects  the  model  while  inside  the 
barrel  during  the  launch  cycle  and  forms  a  gas 
3eal  during  the  burning  of  the  propellant  so  that 
consistent  and  maximum  exit  velocities  are 
achieved.  Upon  exit  from  the  barrel,  the  sabot 
pedals  and  pusher  are  stripped  away  and  only  the 
model  is  allowed  to  continue  through  the 
instrumented  section  of  the  test  facility. 


The  free  flight  tests  were  conducted  at  the 
USAF  Aeroballistic  Research  Facility  (ARP),  Eglin 
AFB,  Florida.  The  ARF  is  an  enclosed  concrete 
structure  used  to  examine  the  exterior  ballistics 
of  various  munitions  while  in  unrestrained 
flight.  The  facility  contains  model  measurement 
equipment,  a  launch  room,  a  blast  chamber  (to 
strip  away  the  3abot),  an  instrumentation  control 
room,  and  the  instrumented  range.  The  range 
atmospheric  conditions  are  closely  monitored.  The 
instrumented  section  of  the  ARF  is  207  meters 
long.  It  has  a  13.90  square  meter  cross-section 
for  the  first  69  meters  and  a  23.80  square  meter 
cross-section  for  the  remaining  length.  There  are 
131  locations  available  as  instrumentation  sites 
and  each  has  a  physical  separation  of  1.52  meters. 
Presently,  50  of  those  sites  are  used  to  house 
orthogonal  shadowgraph  stations.  At  each  of  these 
stations  the  maximum  window  (an  imaginary  circle 
in  which  a  projectile  in  flight  will  cast  a  shadow 
on  both  orthogonal  reflective  screens)  is  2.13 
meters  in  diameter.  The  orthogonal  photographs 
obtained  of  the  model's  shadow  (called 
shadowgrams)  are  subsequently  used  to  determine 
the  spatial  position  and  angular  orientations  of 
the  projectile  at  each  of  the  shadowgraph 
stations.  The  position,  orientation,  and  time 
data  are  then  used  to  reconstruct  the  experimental 
trajectory  which  i3  compared  to  the  governing 
theoretical  equations  of  motion  in  order  to 
extract  the  aerodynamics  of  the  test  projectile. 
This  is  accomplished  by  using  the  Aeroballistic 
Research  Facility  Data  Analysis  System  (ARFDAS).^ 
ARFDAS  is  a  user  friendly  data  reduction  code 
which  enables  the  engineer  to  analyze,  plot,  and 
tabulate  the  aerodynamic  results  from  each  test 
flight.  Figure  11  shows  a  typical  flowfleld 
photograph  taken  from  a  direct  spark  shadowgraph . 


Fig.  3  Photograph  of  model/sabot  packages. 


Fig.  4  Typical  flowfleld  shadowgram  of  a 
small  fin  model  in  flight. 

Prediction  Techniques 


Preceding  the  te3t  flights,  statio  stability 
predictions  were  made  using  both  empirical  and 
computational  methods.  The  empirical  predictions 
were  made  using  the  Projectile  Design  Analysis 
System  (PR0DAS)  .  PR0DAS  is  an  interactive 
program  which  allows  the  engineer  to  create  a 
projectile  by  specifying  it's  material  and 
physical  dimensions.  It  calculates  the  model 
physical  properties  and  compares  that  design  to 
similar  configurations  in  an  extensive  data  base. 
The  program  then  interpolates  the  data  base  to 
provide  estimated  aerodynamic  coefficients  for  the 
created  configuration.  The  following  data  figures 
show  these  aerodynamic  predictions  for  both  the 
small  and  large  delta  fin  configurations 
previously  described. 


□  □ 


The  computational  predictions  were  obtained 
by  using  two  methods.  First,  an  inviscid  space 
marching  Euler  code  was  used  to  predict  the 
stability  characteristics  for  both  configurations 
of  interest  between  Mach  3.0  and  7.0.  The  Zonal 
Euler  Solver  (ZEUS)’  code,  is  a  finite  volume, 
flow  solver  that  uses  shock  fitting  for  the  bow 
shock.  It  was  selected  because  of  its  efficiency 
and  robustness.  The  ZEUS  code  is  especially 
suitable  for  computations  in  which  the  flow  field 
is  completely  supersonic  about  the  body  and 
assumes  that  the  fins  are  infinitely  thin.  The 
code  applies  a  2onal  strategy  which  divides  the 
computational  domain  into  zones  which  can  be 
mapped  separately.  The  zone  edges  are  defined  by 
the  user  to  coincide  with  the  body  and  bow  shock 
surfaces.  The  calculation  of  the  flowfleld  is 
initiated  near  the  nose  tip  of  the  body  and 
marched  aft  in  a  planar  fashion  along  the 
longitudinal  axis  of  the  body  to  its  base.  The 
initial  flowfield  wa3  generated  by  a  separate 
conical-flow  starting  solution.  Three 
computational  zones  were  established  due  to 
symmetry  ccnditior.3.  Two  fin3  were  defined  in  the 
half  plane  as  presented  in  the  sketch  below. 


The  two  fins  are  located  at  the  45°  and  135° 
positions,  respectively,  where  0°  is  defined  in 
the  sketch.  A  9  X  36  grid  was  employed  in  zone  1, 
zone  2  employed  a  18  X  36  grid  and  zone  3  used  a  9 
X  36  grid.  The  grids  in  the  circumferential 
direction  were  equally  spaced,  while  the  36  points 
U3ed  in  the  radial  direction  were  clustered 
between  the  body  surface  and  the  bow  shock  wave. 
These  cells  were  clustered  so  that  grid  points 
were  concentrated  close  to  the  body.  The  step 
size  in  the  marching  direction,  along  t*e  booy 
axis,  was  automatically  varied  as  a  fur  to.,  of 
goemetry,  flow  condition  and  stability  cr  iter's. 

The  Eglln  Arbitrary  Geometry,  Implici  er 
(EAGLE)  Code,  was  used  as  another  analytics 
method  to  predict  static  stability.  This  program 
is  a  multiblock  grid  generation  and  steady  state 
flow  solver  system.  It  combines  a  boundary 
conforming  surface  generation  scheme  ’’ ° ,  a 
composite  block  structure  grid  generation  scheme 
,  and  a  multiblock  implicit  Euler  flow  solver 
algorithm.  0  The  EAGLE  code  is  used  to  obtain 
predictions  for  the  flight  regimes  in  which  the 
flowfield  is  not  necessarily  entirely  supersonic. 
The  computational  grid  developed  for  the  fin- 
stabilized  missile  discussed  herein  consisted  of 
twelve  (12)  three-dimensional  blocks  for  a  total 
of  170,240  grid  point  locations.  Three  regions 
were  defined  in  the  grid  generation  process  as 
existing  between  each  of  the  four  (4)  fins  at  the 
rear  of  the  missile.  Each  of  the  twelve  blocks 
are  defined  by  rotating  the  three  regions  between 
any  two  fins,  a  rotation  of  90  degrees. 


The  grid  points  were  clustered  near  the  body 
surface  using  a  hyperbolic  tangent  distribution 
starting  with  an  initial  spacing  of  .001  in.  Two 
separate  grid  descriptions  were  generated,  one  for 
transonic  flow  (M  =  .8  to  M  =  1.2)  and  another  for 
supersonic  flow  (M  a  1.5  to  M  =  4.5).  The 
dimensions  of  each  grid  field  has  been  optimized 
for  the  characteristics  of  each  flow  regime.  The 
grid  for  the  transonic  regime  (M  x  0.8  to  1.2)  has 
been  dimensioned  to  capture  a  Mach  wave  in  -  sin* 
1/M)  expected  to  propagate  from  a  point  source 
located  at  the  nose-tip  at  Mach  x  1.2  and 
intersecting  a  point  located  on  the  flow-exit 
plane.  The  Mach  wave  will  intersect  the  flow-exit 
surface  at  the  60  percent  point  of  the  exit  plane 
length,  leaving  space  for  the  shock  wave  to  be 
captured.  The  same  reasoning  was  applied  to  the 
supersonic  regime  from  Mach  x  1.5  to  Mach  x  4.5  by 
selecting  Mach  =  2.5  as  the  baseline  condition. 
Figure  5a  Is  the  cross-sectional  view  of  the  grid 
used  for  the  supersonic  series  of  Mach  numbers  (M 
=  1.5  to  4.5)  and  part  b  of  Figure  5  is  a  cross- 
sectional  view  of  the  grid  used  for  the  transonic 
series  of  Mach  numbers  (M  x  .8  to  1.2).  It  was 
found  that  maximum  grid  density  could  be  achieved 
by  using  the  Mach  wave  to  define  the  approximate 
location  of  the  bow  3hock  and  adjusting  the  outer 
boundary  to  allow  the  shock  wave  to  be  captured 
during  the  time  dependent  solution. 


a.  Supersonic  Grid 


b.  Transonio  Grid 


Fig.  5  EAGLE  code  grid. 
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The  aerodynamic  derivatives,  Cm^  and  , 
were  computed  by  differencing  the  normal  force  and 
moment  coefficients  at  zero  and  0.5  deg.  angle  of 
attack.  The  moment  coefficient  was  obtained 
directly  from  the  EAGLE  code  and  the  normal  force 
coefficient  was  computed  using  the  lift  and  drag 
coefficients  as  obtained  from  EAGLE,  (ie,  CN  = 

CL  cos®<  +  CD  sino<).  Since  the  fins  were 
considered  to  be  infinitely  thin  they  were  modeled 
by  using  the  proper  flow-field  boundary  conditions 
between  the  inter-connecting  blocks  in  the  EAGLE 
flow  solver. 

The  results  from  both  the  ZEUS  and  EAGLE  code 
predictions  are  plotted  along  with  the  PRODAS 
predictions  in  the  following  data  figures.  These 
predictions  formed  the  baseline  for  the 
aerodynamic  tests.  The  remainder  of  this  paper 
will  discuss  the  free  flight  test  analysis, 
comparisons  with  the  predictions,  and  the 
resulting  conclusions. 

Experimental/Computational  Results 

Tables  1  and  2  summarize  the  experimental 
results  for  the  small  and  large  fin  model 
configurations  and  include  both  single  and 
multiple  fit  reductions.  The  multiple  fit  is  a 
technique  where  flights  of  the  same  configuration 
and  at  similar  Mach  number  conditions  are 
simultaneously  analyzed.  Thi3  provides  a  common 
set  of  aerodynamics  that  match  each  of  the 
separately  measured  position-attitude-time 
profiles.  This  approach  provides  a  broader  range 
of  angle  of  attack  and  roll  orientation 
combinations  than  would  be  available  from  each 
trajectory  considered  separately.  This  increases 
the  probability  that  the  determined  coefficients 
define  the  configuration's  aerodynamics  over  a 
larger  range  of  possible  trajectories.  There  are 
blanks  in  the  summary  tables  where  that  particular 
value  was  held  fixed  (for  reasons  that  will  be 
explained  in  subsequent  discussions).  This  was 
done  so  that  only  experimentally  determined 
coefficients  would  be  shown  in  the  Tables  1  and  2. 

Pitching  Moment  Derivative 

Figure  6  show3  the  experimentally  determined 
pitching  moment  derivative  for  both  the  small  ar.d 
large  fin  configurations,  respectively.  Note  the 
two  distinct  levels  resulting  from  having  two 
different  center  of  gravity  locations  (as 
described  earlier)  in  each  of  the  plots.  As 
expected,  the  forward  center  of  gravity 
configurations  exhibit  greater  stability  which  is 
exemplified  by  the  more  negative  (parallel  shift) 
derivative  values.  Based  on  this  plot,  one  would 
expect  the  small  fin  (Xcg/1=.52)  model  to  be 
statically  unstable  at  approximately  Mach  s  2.2. 

As  can  be  3een  in  Figure  7,  which  plots  the 
angular  motion  versus  travel  for  a  model  fired  at 
this  very  condition,  this  is  precisely  what 
occurred.  This  particular  flight  was  the  only  one 
purposely  launched  at  such  a  condition  because  we 
did  not  want  to  risk  damage  to  the  ARF 
instrumentation  systems. 

The  majority  of  the  flights  experienced 
relatively  small  angular  motions,  therefore  the 
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Fig.  6  Pitching  moment  derivatives. 
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Fig.  7  Angular  motion  plot  of  statically 
unstable  flight. 


determination  of  non-linear  pitching  moment 
coefficients  was  not  attempted  for  these. 

However,  for  flights  that  did  have  a  significant 
angle  of  attack,  an  attempt  to  fit  the  non-linear 
terms  was  tried  but  these  attempts  were  not 
successful.  The  small  angular  motions  for  many  of 
these  flights  also  made  it  more  difficult  to 
obtain  accurate  zero  angle  of  attack  coefficients, 
such  as  the  first  order  pitching  moment 
coefficient.  Note  that  some  Cm^  values  in  the 
tables  have  been  deleted  for  this  very  reason. 
Additionally,  notice  that  the  Mach  number 
dependent  pitching  moment  term,  Cm^,  can  be 
easily  obtained  from  the  slopes  of  the  curves  in 
Figure  6.  These  values  are  approximately  40  and 
30  for  the  small  and  large  fin  models 
respectively.  It  should  be  noted  that  the  angular 
motion  fit  (to  the  flight  trajectories)  wa3  not 
very  dependent  on  having  an  accurate  value  for 
thi3  coefficient. 


Figure  6a  displays  the  results  for  Cm  as 
computed  by  ZEUS  and  PRODAS  for  the  small ''fin 
model.  The  ZEUS  and  PRODAS  results  were  computed 
for  the  X0  /I  i  .43  location  and  not  for  the  Xcg/1 
*  .53.  The  PRODAS  results  are  clearly 
underpredicting  the  C  data  for  the  free-flight 
measurements  and  the  ZEUS  code  appears  to  predict 
these  results  better  if  the  computations  could  be 
extended  to  the  lower  Mach  numbers  Indicated  by 
the  test  data.  However,  since  ZEUS  is  a 
supersonic  marching  code  and  the  available  data  is 
in  the  transonic  range.  Mi  8  to  M  i  1.84,  the 
ZEUS  code  would  become  numerically  unstable 
because  of  the  hyperbolic  nature  of  the  supersonic 
Euler  equations.  Therefore,  direct  computation  in 
this  range  was  not  attempted  using  the  ZEUS  code. 
The  apparant  discrepancy  between  the  free-flight 
data  and  the  theoretical  results  indicate  that  fin 
size  may  have  been  a  major  factor  in  the 
difference.  It  is  indicated  that  the  boundary 
layer  thickness  at  the  end  of  the  missile  is  thick 
enough  to  render  the  Prodas  solution  ineffective 
for  the  3mall  fin  model.  Not  much  can  be  3aid 
about  the  Zeus  results  because  the  predictions  do 
not  extend  far  enough  into  the  region  covered  by 
the  free-flight  data. 


Figure  6b  shows  the  results  for  C_  as 

UJ  cf 

predicted  by  EAGLE,  ZEUS  and  PRODAS  for  the  large 
fin  configuration.  The  predictive  techniques  were 
executed  for  a  Xcg/1  reference  location  equal  to 
0.43.  The  EAGLE  code  and  the  PRODAS  results  are 
reasonably  close  over  the  entire  range  from  M  i  .8 
to  M  i  4.5  as  displayed  on  the  plot.  It  is 
significant  that  such  good  correlation  with  the 
free-flight  data  can  be  attained  over  a  range  of 
Mach  numbers  that  extends  from  Mi. 8  to  Mi  4. 5. 
Free  Flight  data  is  under-predicted  to  a  small 
degree  by  all  of  the  methods,  but  the  trend  in  the 
data  is  captured  by  all  of  the  CFD  techniques.  It 
is  interesting  to  note  that  the  ZEUS  code  seems  to 
predict  C_  better  than  either  of  the  other  two 

lUO( 

methods  at  the  higher  Mach  numbers.  The  ZEUS  data 
however,  converges  to  EAGLE  and  PRODAS  results  as 
the  Mach  number  is  reduced  below  M  i  3. 


Normal  Force  Derivative 


Because  the  force  coefficient  effects  are 
more  difficult  to  extract  than  moment  derivatives 
from  flights  with  little  angular  motion,  the 
normal  force  coefficient  was  understandably  harder 
to  obtain  from  these  flights.  However,  this  was 
not  unexpected,  and  was  the  reason  that  the  two 
different  center  of  gravity  model  configurations 
were  fabricated.  Figure  8  shows  the  determined 
normal  force  coefficient  derivative  using  the 
multiple  fit  method  previously  discussed  plotted 
versus  Mach  number  for  the  small  and  large  fin 
configurations.  Also,  shown  in  this  figure,  are 
empirically  determined  values  obtained  using  the 
technique  shown  in  Figure  9.  This  figure  plots 
the  experimentally  determined  pitching  moment 
versus  Xcg/1  for  Mach  number  i  1.2  (  small  fin 
model ) . 
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Fig.  8  Normal  force  derivatives. 
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The  slope  of  the  curve  generated  by  the  two 
center  of  gravity  locations  is  a  measure  of  the 
normal  force  coefficient  derivative  obtained  from 
the  following  relationship: 


CN  * _ * _ 

<*  (1/d)  (Xcg/l-Xcp/1) 

where  Xcp/1  corresponds  to  where  the  curve  crosses 
the  X-axi3  at  Cnc<  =  0.  Using  this  approach  and 
performing  the  calculation  results  in  a  CN(X 
estimate  of  9.06  for  the  small  fin  configuration 
at  both  M=1. 2  and  2.2.  The  large  fin 
configuration  resulted  in  a  CNlJ<  value  of  9.45  at 
H-2.2  and  a  value  of  7.29  at  M=3«35.  These  values 
were  confirmed  further  by  inputing  them  in  the 
analysis  program,  for  each  of  the  flights  where 
the  normal  force  was  not  experimentally  accurately 
determined,  and  holding  that  value  fixed.  In  each 
case,  the  probable  error  of  fit  (in  swerve  and 
angle)  were  not  adversely  affected  and  in  most 
oases  the  fits  benefited  slightly.  In  addition, 
note  that  the  multiple  fit  values  shown  in  Figure 
8  converged  to  solutions  with  enclose  to  these. 
It  is  believed  tha  the  scatter  in  the  fitted 
values  is  caused  by  the  relatively  small  angles  of 
attack  experienced  during  the  flights  resulting  in 
very  small  lateral  movement  of  the  model.  The 
lateral  movement  is  the  source  of  the  fitted 
normal  force  derivatives.  It  shall  be  noted  that 
blank  values  for  C jg  in  the  summary  tables 
indicate  that  the  empirically  determined  value  was 
used  and  held  constant  during  the  fitting  process. 
Because  of  the  relatively  large  data  scatter  all 
that  can  be  said  for  the  comparison  between  the 
experimental  data  and  the  theoretical  predictions 
are  that  the  levels  are  in  general  agreement. 


Center  of  Pressure  Location 

Center  of  pressure  location  a3  a  function  of 
Mach  number  are  presented  in  Figure  10  for  the 
large  and  small  fin  configurations,  respectively. 

The  experimental  data  shown  in  this  figure  wa3 
obtained  from  the  previously  discussed  CBt<and  CN 
results  using  the  following  relationships:  '  * 


Fig.  10  Center  of  Pressure  Location. 
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The  PRODAS  and  EAGLE  prediction  for  the  large 
fin  configuration  compare  well  over  the  entire 
range  of  Mach  numbers.  ZEUS  data  converges  with 
the  EAGLE  and  PRODAS  results  in  the  region  from 
M=3  to  M=4.5,  (see  figure  10b).  The  same  level  of 
convergence  for  ZEUS  and  PRODAS  is  observed  in 
figure  10a  for  the  small  fin  model.  The  free- 
flight  data  for  both  fin  configurations  show  good 
agreement  between  the  theoretical  and  experimental 
results . 

Axial  Force  Coefficient 


The  axial  force  coefficient  (CXQ)  is  one  of 
the  easiest  and  most  accurate  coefficients  to 
obtain  from  free  flight  aeroballistic  tests. 

Figure  11  shows  the  experimentally  obtained  zero 
yaw  axial  force  coefficients  plotted  versus  Mach 
number.  The  EAGLE  and  PRODAS  predictions  are  also 
plotted  on  this  figure.  The  EAGLE  prediction  has 
been  corrected  by  adding  skin  friction  drag  and 
base  drag  to  the  wave  drag  component  computed  by 
the  inviscid  EAGLE  Flow  solver. 
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The  relationship  chosen  to  model  the  turbulent 
skin  friction  coefficient  has  the  following  form: 


CDf 


0.455  /Jw_\ 

(log  Rx)  2-58  \  Sb  J 


(3) 


Where  Sw/Sb  is  the  rafio  of  the  wetted  area  to  the 
base  area  of  the  model'.  Equation  (3)  is  valid  for 
Reynolds  numbers  ranging  from  5  X  10’  to  1  1  10  . 
The  flow  over  the  entire  surface  of  the  model  is 
assumed  to  be  turbulent.11 

The  base  drag  component  is  obtained  using  one 
of  two  empirical  relations  depending  on  whether  or 
not  the  flow  is  supersonic  or  subsonic: 


Subsonic,  CDb  = 

0.029/  CDC 

(4) 

Supersonic,  Cpb 

z  K  (1.43/M2) 

(5) 

Where  K  is  a  function  of  Mach  number.11 


Figure  11b,  which  presents  the  results  of  the 
large  fin  analysis,  shows  good  correlation  of 
experimental  free-flight  data  and  corrected  EAGLE 
results  in  the  middle  of  the  Maoh  number  regime. 
Figure  12  presents  the  various  components  (Cp., 
CDb,  and  C^)  that  make  up  the  total  theoretical 
drag  coefficient  and  Figures  13  and  14  present  the 
density  contours  for  a  typical  transonic  and 
supersonic  EAGLE  code  analysis.  Two  different 
grid  densities  were  used  for  these  computations  as 
previously  discussed,  and  represents  an  attempt  to 
tailor  the  grid  to  the  desired  flow  field. 


Fig.  12  Drag  Components. 


Fig.  11  Axial  force  coefficient. 
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Fig.  13  Transonic  Density  Contour  (<x  =  0.5  deg.) 


Fig.  1 4  Supersonic  Density  Contour  (ot=  0  deg.) 


Roll  Damping  and  Pitch  Damping  Derivative 

The  experimentally  determined  roll  damping 
derivative  (C,  )  and  the  pitch  damping  derivative 
(C  )  are  tabulated  In  Tables  1  and  2.  Similarly, 
as  for  the  other  coefficients,  a  blank  in  the 
summary  tables  for  these  parameters  means  that 
they  were  held  constant  during  the  analysis  of 
that  flight.  These  results  (both  Clp  and  C  -)  are 
not  plotted  herein  because  the  ZEUS  and  EAGLE  code 
predictions  were  not  available.  The  ZEUS  code  is 
unable  to  accomplish  these  predictions  and  the 
predictions  from  the  EAGLE  code  were  not 
accomplished  because  these  predictions  require  an 
unsteady,  time  dependant  calculation  which  at 
present  is  not  included  in  the  routine. 

Upon  inspection  of  the  C^p  data  shown  in 
Tables  1  and  2  these  derivatives  appear  to  have 
considerable  scatter.  This  is  believed  to  be 
caused  by  the  very  low  roll  rates  experienced 
during  the  flights  (i.e.,  the  fins  were  not 
canted).  It  is  significantly  more  difficult  to 
extract  C^p  from  flights  experiencing  low  roll 


rate3  because  there  is  little  rolling  motion  to  be 
damped.  This  is  analogous  to  determining  the 
pitch  damping  derivative  from  flights  with  very 
low  angular  motion.  The  Cog  data  tabulated  in 
Tables  1  and  2  indicate  that  both  the  large  and 
small  fin  models  are  dynamically  stable.  The  vast 
majority  of  the  experimentally  obtained  Cm(J  values 
lie  between  -  2000  and  -  t000  per  radian  with  no 
apparant  effect  due  to  fin  size. 

Conclusions 

A  comparison  of  experimental  free-flight  data 
with  three  different  prediction  routines  have  been 
accomplished  for  a  slender  fin  stabilized  missile 
configuration.  These  comparisons  have  been 
accomplished  in  the  Mach  number  range  of  1.0  to 
5.0  and  in  general  are  in  good  agreement.  All 
three  prediction  routines  indicates  that  the 
static  stability  decreases  rapidly  with  increasing 
Mach  number  and  are  consi3tant  with  the 
experimental  data.  Although  the  present  analysis 
does  not  identify  the  cause  of  this  decreasing 
static  stability,  it  is  apparantly  not  caused  by 
flow  viscous  effects  or  the  loss  of  fin 
effectiveness  due  to  the  boundary  layer  thickness. 
This  is  apparant  because  both  the  ZEUS  and  EAGLE 
codes  are  invlsid  routines  and  they  predict  this 
loss  of  static  stability. 

It  is  also  of  interest  to  note  that  the  ZEUS 
code  predicts  that  the  forward  shift  of  the 
center-of-pressure  ccase3  near  Mach  five 
indicating  that  these  high  fineness  ratio,  tall 
control  configurations  may  have  applications  in 
the  hypersonic  region.  Unfortunately  the  present 
experimental  data  does  not  extend  far  enough  into 
this  region  (above  M=5)  to  confirm  this  trend.  It 
is  hoped  that  this  wiU  be  accomplished  in  future 
tests.  Also,  the  EAGLE  code  corrected  for  skin 
friction  and  base  drag  reasonably  predicted  the 
axial  force  coefficient.  It  is  believed  however, 
that  this  prediction  could  be  further  improved  by 
U3ing  more  sophisticated  correction  techniques. 
This  's  also  an  *r»a  were  future  work  will  be 
accomplished . 

Although  these  tests  were  restricted  to 
researching  the  stability  characteristics  of  delta 
fin  configurations,  the  test  program  is  continuing 
to  include  other  parameters.  For  example,  the 
next  phase  (to  be  accomplished  during  the  summer 
of  89)  will  investigate  the  potential  advantage 
associated  with  moving  mere  of  the  fin  area 
outside  the  boundary  layer.  A  clipped  delta  fin 
configuration  has  been  fabricated  which  has  the 
same  fin  area  as  the  large  delta  fin  models,  while 
maintaining  the  fin  span  of  the  small  delta  fin 
models  (described  in  this  paper).  In  addition,  an 
inverted  fin  model  (where  the  root  chord  is  very 
small  compared  to  the  tip  chord)  has  been 
constructed  which  also  has  the  same  area  and  span 
constraints  described  above.  It  is  hoped  that  the 
static  stability  of  these  configurations  will  be 
improved  by  moving  more  of  the  fin  outside  the 
boundary  layer  thus  improving  the  model  fin 
effectiveness.  In  addition  to  these  near  term 
tests,  it  is  expected  that  testing  will  continue 
into  future  years  where  additional  parameters  will 
be  Investigated,  such  as  fin  longitudinal 
location,  nose  shape  variations,  and  even  flare 
variations . 
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Table  1  Experimental  results  (large  fin) 
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-*2.68 

-3182 

-5-93 

0.*3 

1.828 

*.9 

0.668 

11.2* 

-40.86 

-i860 

-3.09 

0.43 

iHCttfci  rlfi 


1.80* 

TT 

0.725 

11.77 

-53.U 

-269* 

0.43 

1.856 

2.* 

0.639 

7.30 

-23-86 

-25*7 

-2.89 

0.52 

1.8*0 

1.1 

0.655 

9.15 

-*2.22 

-2822 

-6.  *1 

0.*3 

11.  Hoerner,  Sighard,  "Fluid  Dynamic  Drag", 
Hoerner  Fluid  Dynamics,  1 9 65 . 
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